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■ Abstract 

<N 

^ | We consider for j = ~, |, . . . a spherically symmetric, static system of (2j + 1) 

. Dirac particles, each having total angular momentum j. The Dirac particles interact 

pLn ' via a classical gravitational and electromagnetic field. 

The Einstein-Dirac-Maxwell equations for this system are derived. It is shown 
that, under weak regularity conditions on the form of the horizon, the only black hole 
solutions of the EDM equations are the Reissner-Nordstrom solutions. In other words, 
the spinors must vanish identically. Applied to the gravitational collapse of a "cloud" 
of spin-i-particles to a black hole, our result indicates that the Dirac particles must 
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eventually disappear inside the event horizon. 



00 ! 1 Introduction 

ON 

In [TJ [2] , particle- like solutions of the Einstein-Dirac-Maxwell (EDM) equations were con- 
structed for a static, spherically symmetric singlet system. It was found that the solutions 
in a given state (i.e. in the ground state or in any fixed excited state) cease to exist if the 
rest mass m of the fermions becomes larger than a certain threshold value m s . The most 
natural physical interpretation of this observation is that for m > m s , the gravitational in- 
teraction becomes so strong that a black hole would form. This suggests that there should 
be black hole solutions of the coupled EDM equations for large fermion masses. The work 
[3], however, indicates that this intuitive picture of black hole formation is wrong. Namely, 
it was proved in [3] that the Dirac equation has no time-periodic solutions in a Reissner- 
Nordstrom black hole background, even if the Dirac particles have angular momentum 
and can thus, in the classical picture, "rotate around" the black hole. This implies that 
if we neglect the influence of the Dirac particles on the gravitational and electromagnetic 
field, there are no black hole solutions of the EDM equations. 

In order to understand if and how Dirac particles can form black holes, we study in 
this paper the fully coupled EDM equations. We do not assume the Dirac particles to be 
in a spherically symmetric state; they are allowed to have angular momentum j. However, 
we arrange 2j + 1 of these particles in such a way that the total system is static and 
spherically symmetric. In the language of atomic physics, we consider the completely 
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filled shell of states with angular momentum j. Classically, one can think of this multiple- 
particle system as of several Dirac particles rotating around a common center such that 
their angular momentum adds up to zero. Since the system of fermions is spherically 
symmetric, we get a consistent set of equations if we also assume spherical symmetry for 
the gravitational and electromagnetic field. This allows us to separate out the angular 
dependence and reduce the problem to the analysis of a system of nonlinear ODEs. 

We prove analytically that, under weak regularity conditions on the form of the horizon, 
the black hole solutions of our coupled EDM equations are either the Reissner-Nordstrom 
solutions (in which case the Dirac wave functions are identically zero) , or the event horizon 
has the form of the extreme Reissner-Nordstrom metric. In the latter case, we show 
numerically that the Dirac wave functions cannot be normalized. Thus our Einstein- 
Dirac-Maxwell system does not admit black hole solutions. Our results show that the 
study of black holes in the presence of Dirac spinors leads to unexpected physical effects. 
Applied to the gravitational collapse of a "cloud" of Dirac particles, this is a further 
indication that if an event horizon forms, the Dirac particles must eventually disappear 
inside this horizon. 

The methods used in this paper are quite different from those in [3]. Namely, in 
contrast to [3], we do not derive "matching conditions" for the spinors across the horizon. 
We work here only with the equations outside the event horizon, and the proof relies on 
the nonlinear coupling of the spinors into the Einstein-Maxwell equations. 



2 The Spherically Symmetric Multi-Particle System 

The gravitational field is described by the static, spherically symmetric Lorentzian metric 
in polar coordinates 

ds 2 = g ij dx i dx j = dt 2 - -dr 2 - r 2 d$ 2 - r 2 sin 2 dip 2 (2.1) 

with positive functions A = A{r) and T = T(r). We consider as our space-time the region 
r > p > outside a ball of radius p around the origin. The physical situation we have in 
mind is that the surface r = p is the event horizon of a black hole. We assume the metric 
to be asymptotically Minkowskian, 

lim Air) = 1 = lim T(r) . (2.2) 

r— >oo r— »oo 

The electromagnetic field is described by a potential A of the form A = (—0,0), where <p 
is the Coulomb potential <p = (j)(r). 

In direct generalization of the situation in the Reissner-Nordstrom background [3] , the 
Dirac operator G takes the form 




where 7*, j r , 7 , and 7^ denote the 7-matrices of Minkowski space in polar coordinates, 

l l = 7° 

7 r = 7 1 cos 1? + j 2 sint? cos <p + 7 3 sin?? sin Lp 
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7^ = - ( — 7 1 sini? + 7 2 cost? cos cp + 7 3 cos 1 !? sin^j 

7 V = (—j 2 sin 99 + 7 3 cos(/?) 

rsmv V / 

As with the central force problem in Minkowski space [3|, this Dirac operator commutes 
with: a) the time translation operator idt, b) the total angular momentum operator J 2 = 
(L + S) 2 , c) the z-component of total angular momentum J z , and d) with the operator 
7° P (where P is parity). Since these operators also commute with each other, we can 
write any solution of the Dirac equation as a linear combination of solutions which are 
simultaneous eigenstates of these operators. We denote this "eigenvector basis" for the 
solutions by 

1 3 

^%uj wi th c = ±, j = -,-,..., k = -j,-j + wGlR; (2.4) 

the eigenvalues are 

J 2 * c jku 
J* 

7°P $± 

For the functions ^j ku , the Dirac equation reduces to ordinary differential equations in the 
radial variable r. The quantum number k describes the orientation of the wave function 
^jkuj i n s P ace ) an d thus spherical symmetry of the Dirac operator implies that the radial 
Dirac equation is the same for all values of A;. In order to build up our many-particle 
system, we take, for given c, j, and uj, one solution of the radial Dirac equation and 
consider the system of the 2j + 1 particles 

^ ' jkui ) ^ 3i ■ ■ ■ i j 

corresponding to this radial solution. Using the formalism of many-particle quantum 
mechanics, we can describe the fermions with the Hartree-Fock state 

if HF = . A^ c i. ■ 1 1 A • • • A ^ c i. (2 5) 

j k=~j u> ' x ^ j k=—j+l u> ' x 3 k=j u> I 

For simplicity, we usually avoid this formalism here; it is easier to just work with an 
orthonormal basis of the one-particle states. 

For clarity, we point out that each fermion has non-zero angular momentum and is 
thus not in a spherically symmetric state. Nevertheless, the system of 2j + 1 particles is 
spherically symmetric; this can be verified in detail as follows: The Hartree-Fock state 
(12. 5ft is an eigenstate of J z ; namely 

(Jz *i k=-j J A ^ k= _ j+1 w A • • • A ^ k=j u 

k=~j u A (J z ^ fc= _ j+1 J A • • • A ^ k=j u 
+ ... + ^ k= _ 3 u A k= _ j+1 u A • • • A ( J z V? k=jw ) 

]T fc* HF = . 
k=-j 



U Vjkuj 
j(j + 1) Vjkuj 
jkuj 



k tyf 



_,_ i 1 for j + i even 



' ^ jkw X 1 -1 for j + i odd 
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Similarly, we can apply the "ladder operators" J± = J x ± iJ y to the Hartree-Fock state, 

J± ^ HF = ( J± 4* k= _ 3 J A ^ c k= _ j+l u A • • • A ^ k=j u 

k=-j cu A ( J± * • k=-j+lJ A • • • A *i fc=j - 
+ • • • + ^ *=_,- w A ^ c fe= _ J+1 w A • • • A ( J± k=j J . (2.6) 

After substituting the relations 

J±*jk» = ^Jj(j + l)-k(k±l) Vj k±lLU , 

the anti-symmetry of the wedge product yields that each summand in (|2.6p vanishes. We 
conclude that 

J* HF = . 

Since the total angular momentum operator J is the infinitesimal generator of rotations, 
this implies that ^ HF is spherically symmetric. 

For a physically meaningful solution of the Dirac equation, the wave function must be 
normalized. The normalization integral for the wave functions ^f^u over the hypersurface 
TL = {t = const, r > p} is 

where v is the future-directed normal of TC, and where d/1-n denotes the invariant measure 
on TL induced by the Lorentzian metric. For a normalized solution of the Dirac equation, 
(12 .7j) gives the probability for the particle to be in the region r > p outside the ball of 
radius p centered at the origin. It seems tempting to demand that, for a normalizable 
solution of the Dirac equation, the integral (|2.7j) must be finite. However, as explained in 
[3], the normalization integral inside the event horizon (which we do not consider here), 
is not necessarily positive, and it might happen that an infinite contribution near r = p 
in (|2.7p is compensated by an infinite negative contribution inside the horizon. Therefore 
we only demand that the normalization integral away from the horizon is finite; namely 



^jku^^jku Vj d ^n < oo for every r > p. (2.8) 

{t=const, r>ro} 

We remark that the singlet state of [lj [2] can be recovered from our multi-particle 
system by considering the case j = i. 

3 The EDM Equations 

We now derive the system of differential equations. We begin by separating out the angular 
and time dependence in the Dirac equation. We choose the wave functions VE^^ of the 
previous section in the explicit form 
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where are two-component radial functions, and where X + 1> 3 = k-> !>•••> & 

—j, + denote the 2-spinors 



k j + l-k k-i I l 

K i+h V 2j + 2 j+5 I 




(Y ; m are the spherical harmonics). The functions X-+± are eigenvectors of the operator 
K = aL+ 1, 

,fc -W„- , ^ 



Using the relations between the angular momentum operators (see |3j for details), the Dirac 
equation {G — m)^j kuJ with G as in (|2.3|) reduces to the ordinary differential equation 




(w-e4>)T\ ~ nM± ^l JO 




• (3-3) 



Since this equation is independent of A;, we will in the following omit the index k and 
simply write The normalization integral (|2.8p for the wave functions takes the form 



f°° \J x 

j \^ C ju){ r )\ 2 dr < oo for every tq > p. (3.4) 
Jr A 

The Dirac equation (13.3H implies that the "radial flux" function (see [3J) 



m ■■= $u r ) [i o) ^ (r) (3 - 5) 

is a constant, as is verified by a short computation. Since I^^J 2 > and since the metric 
is asymptotically flat, the normalization integral (|2.7p can be finite only if this constant 
is zero. We can thus assume that F vanishes identically. This means that the product 
^juii^jui2 (constructed from the two components of is a real function. For a given 
radius r$ £ (p, oo), we can thus arrange with a constant phase transformation 

- e ia ^{r) , a £ IR, 

that both ^^1(^0) an d ^jaj2( r o) are real. Since all the coefficients in the Dirac equation 
(|3,3p are real, it follows that the spinors <& C j u {r) are real even for all r £ (p, 00). We denote 
these two real components of ^ u by a and (3. 

Next we must calculate the total current and energy-momentum tensor of the Dirac 
particles. With our explicit formulas (|3.ip and f)3.2|) for the angular dependence of the 
wave functions, the anti-symmetrization in the Hartree-Fock state (j2.5j) is trivial. One 
obtains that the electromagnetic current of the multi-particle system is simply the sum of 
the currents of all states ^j kuj , 

j 



i k = * HF G k * HF = £ *° jku G k * c jku} 

k=-j 
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Because of spherical symmetry, the angular components j and of the Dirac current 
vanish. The "sum rule" 

k=-j 

yields, after a straightforward computation, that 

k=-j 

To calculate the radial current j r , we will need the formula [3j Eqn. (3.7)] 

where a r is the Pauli matrix in the radial direction, 

a r = a 1 cos $ + a 2 sini9 cos^ + <r 3 sin-!? sin 93 . 
Using (|3.7p . (13. 6p . and the fact that the radial flux F(r) vanishes, we obtain 

f = E ^^(x)^ = . 

Similar to the total current, the energy-momentum tensor of the multi-particle system is 
simply the sum of the energy- momentum of all states ^j^', it has the general form 

1 t 

T ab = -Re E Vjku, ( iG - (°b ~ ie ^b) + iG b (d a - ieA a )) * c jhl) (3.8) 

k=-j 

(this formula is obtained by varying the classical Dirac action; see e.g. PQ). The following 
calculation depends on whether the index c of the wave functions in (|3.8p is c = + or 
c = — . We use the ±/=F-notation, whereby the upper and lower choice correspond to the 
cases c = + and c = — , respectively. From spherical symmetry, Tj, T^, T£, T^, and 
must vanish. An easy computation using the sum rule (|3.6[) and the Dirac equation (|3.3p 
gives 

_ (^)i! (a2+ , 2) ^+i 

r z 47r 

r z 47T r J 47r r z 47r 

The calculations of T$ and T£ are slightly more difficult. First, spherical symmetry yields 
that 

T$ = T£ = ^Re E + ^ . 

k=-j 

The formula 
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allows us to rewrite the angular derivatives using the angular momentum operator L. This 
gives 

k=-j 



J33 a/3 T 



j 



. Re V x fc ^i ^X fc j_i - X fc i ^X fc i 

fc=-j 



We now use the fact that the spinors X <i are eigenvectors of the operator BL = K — 1, 

^ ± 2 



and carry out the /c-summation with (|3.6p to get 

Finally, we put the obtained formulas for the Dirac current and energy-momentum ten- 
sor together with the Maxwell energy-momentum tensor [2] into the Einstein and Maxwell 
equations 

R*-Ir5} = SttTj , V t F kl = 4vre/ . 

The Einstein equations reduce to two first-order equations for A and T; similarly, Maxwell's 
equations simplify to one second-order equation, and we end up with the following system 
of EDM equations: 



2j + l 

a. — yytju — txp 

2j + l 
2r 

r A' = l-A- 2(2j + l)(w - e(f>)T 2 (a 2 + f3 2 ) - r 2 'AT 2 \4>'\ 2 (3.11) 

r 

T 



Via' = ±^—— a - ({oj - e<p)T + m) (3 (3.9) 
2r 

v 7 !/?' = (fa> — ed>)T — m)a^ -^r— P (3.10) 

2r 

uj>)i 

T' (2i + I) 2 

It A — = A - 1 - 2(2 j + l)(a; - ecfyT 2 (a 2 + /? 2 ) ± 2 W 7 T a/? 



+2(2j + 1) mT (a 2 - /T) + r AT \fi\ z (3.12) 
r 2 A0 " = _( 2i + i) e ( a 2 +/3 2 ) _ ^rA + ^A^ + ^A^' . (3.13) 

The two cases for the signs ±/=p correspond to the two values c = ± for the fermionic wave 
functions ^jkuj- Equations (|3.9p and (|3.10p are the Dirac equations (|3.3p . The Einstein 
equations are (|3.11|) and (|3.12|) . and (|3.13|) is Maxwell's equation. According to (|3.4|) . the 
normalization condition is 

(a 2 + /3 2 )- J 

'ro 



/ (a 2 + /3 2 ) — — dr < oo for every r > p. ( 3 -14) 

Jrn A 



We remark that the system f|3.9[) — fj3.13p . (|3.14p has particle-like solutions, which can 
be constructed numerically using the methods in pp. The mass-energy spectrum of the 
solutions has the same qualitative behavior as for the two-particle EDM system [2J. 
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4 Non-Existence Results 



We want to investigate black hole solutions of the system ()3.9j) — ()3.13j) . This means, more 
precisely, that we consider solutions of f)3.9j) (|3. 13|) defined outside the ball of radius 
p around the origin which are asymptotically flat, (|2,2p . and satisfy the normalization 
condition (|3.4p . We assume that the surface r = p is an event horizon; i.e. the function 
A(r) tends to zero for r \ p, whereas T(r) goes to infinity in this limit. We make the 
following assumptions on the form of the horizon: 



(I) The volume element y|det gij\ = r 2 A 2 T 1 is smooth and non-zero on the horizon, 

i.e. 

T- 2 A-\ T 2 A € C°°([p,oo)) . 
This assumption is sometimes called: the horizon is regular. 

(II) The strength of the electromagnetic field is given by the scalar F^ i™ = —2\<p'\ 2 AT 2 

with the electromagnetic field tensor Fij = diAj — djAi. We assume this scalar to 
be bounded near the horizon; thus in view of (I) we assume that 

\4>'{r)\ <ci , p<r < p + e 1 

for some positive constants c±, e\. 

(III) The function A(r) obeys a power law, i.e. 

A{r) = c(r-p) s + 0((r-p) s+1 ) , r>p (4.1) 

with positive constants c and s. 

If assumptions (I) or (II) were violated, an observer freely falling into the black hole would 
feel strong forces when crossing the horizon. Assumption (III) is a technical condition 
which seems general enough to include all physically relevant horizons. For example, the 
Schwarzschild horizon has s = 1, whereas s = 2 corresponds to the extreme Reissner- 
Nordstrom horizon. However, assumption (III) does not seem to be essential for the state- 
ment of our non-existence results; with more mathematical effort, it could be weakened 
or perhaps even omitted completely. We now state our main result: 

Theorem 4.1 The black hole solutions of the EDM system 13.9\) - [3.13\) satisfying the 
regularity conditions (I), (II), and (III) either coincide with a non- extreme Reissner- 
Nordstrom solution with a = = (3, or s = 2 and the following asymptotic expansions 
hold near r = p: 



A(r 
T(r 
cj)(r 

a(r 
(5{r 



A (r - p) 2 + 0((r - pf) (4.2) 

ro(r-p)- 1 + 0((r-pf) (4.3) 

- + <h{r-p) + 0((r-p) 2 ) (4.4) 
e 

a (r - pf + 0{{r - p) K+l ) (4.5) 

Po(r- P r + 0((r- P r +1 ) (4.6) 



with positive constants Aq, Tq and real parameters 4>o, ao, (3q. The power k must satisfy 
the constraint 
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and the spinor coefficients ag and (3q are related by 



a [ JTqk ± ) = -/3 (m - e^To) , (4. 



2p 

where '± ' refers to the two choices of the signs in Ii3.9\) - (3.13\) . 

We begin the analysis with the case that the power s in (|4.1|) is less than 2. 

Lemma 4.2 Assume that s < 2 and i/iai (a, /3,A,T,(f)) is a black-hole solution where the 
spinors (a, (3) are not identically zero. Then the function (a 2 + j3 2 ) is bounded both from 
above and from below near the horizon, i.e. there are constants c, e > with 

c < air) 2 + I3{rf <-, p <r < p + e . (4.9) 

c 

Proof: The Dirac equations (|3.9|) . (|3.1U|) imply that 




< (4m 2 + (2j r + 2 1)2 V(« 2 + /? 2 ) ■ (4.10) 



Since (a, j3) is a non-trivial solution, the uniqueness theorem for solutions of ODEs implies 
that (a 2 + /3 2 )(r) is non-zero for all p < r < p + e. Thus we can divide Eq. (|4.1U|) by 
\[A~ (a 2 + (3 2 ) and integrate. This yields the bound 



\og((a z + P 2 )(p + e)) - log((^ + (5 2 ){r)) 

/7'+- : / „ (2/ + li 2 Y 



^A-tM^'+W + lgj^ (4u; 



Since s < 2, we see that the integrand in (14.111) is integrable aX r = p. Thus the right-hand 
side of (|4.1ip is majorized by 



P+£ A-Ht)(,rn 2 + ^^-Xdt 



t2 

so we may take the limit r \ p in (|4.1ip to obtain the estimate (|4.9p . 



Proposition 4.3 For < s < 2, the only black hole solutions of the system i3.9\) ~ [3A3\) 
are the non- extreme Reissner-Nordstrom solutions. 

Proof: We shall assume that we are given a black hole solution which is not the Reissner- 
Nordstrom solution, and obtain a contradiction. In this case, the spinors (a, (3) are not 
identically zero, so we may apply Lemma 14.21 and conclude that the spinors are bounded 
near r = p. 
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We first consider the differential equation for AT 2 . The Einstein equations (|3.1ip and 
dSHSD give 

r^{AT 2 ) = -4(2j + 1) (u - e<P) T 4 (a 2 + (3 2 ) ± 2 (2j + 1)2 T 3 a/3 
ar r 

+2(2j + l)mT 3 (a 2 -/? 2 ) . (4.12) 

According to the regularity condition (I), the left and thus also the right side of this 
equation is smooth. Since the spinors are bounded away from zero near r = p, and since 
T — > oo as r \ p, we see that 

lim (uj-e6(r)) = . (4.13) 

p<r^p 

We write Maxwell's equation f)3. 13|) in the form 

<t>" = A (2j+ 9 1)e (a 2 +/3 2 ) - _-L-[r 2 yZT]V . (4.14) 



According to the regularity condition (I) , the square bracket in (|4.14|) , and thus the whole 
coefficient of </>', is a smooth function. However, the factor A -1 in the first summand in 
(|4.14|) blows up on the horizon. If s > 1, the singularity of A^ 1 is not integrable. This 
implies that \(p'\ is unbounded on the horizon, contradicting the regularity condition (II). 
We conclude that s < 1. We can then integrate Eq. (|4.14|) and obtain the local expansion 
around the horizon 

<j/{r) = Cl (r-p)- s+1 + c 2 + 0{(r-p)- s+2 ) , 

where c 2 is an integration constant. Integrating once again and using (|4.13p yields the 
expansion 

<j,(r) = Cl (r - p)- s+2 + c 2 {r-p) + - + 0((r - p)' s+3 ) . (4.15) 

e 

Finally, we substitute the expansion (|4,15p into the A-equation (|3.1ip . Since the 
functions A and r 2 AT 2 \<p'\ 2 are bounded near the horizon, and since (|4.15p implies that 
(u -e<fr) = 0(r - p), whereas T 2 (a 2 + (3 2 ) ~ (r - p)~ s with s < 1, the right side of (|3.1ip 
is bounded in the limit r \ p. However, the left side diverges like rA'(r) ~ (r — pY^ 1 . 
This is a contradiction. ■ 

It remains to consider the case s > 2. 
Lemma 4.4 If s >2, 

lim (r - p)~i (a 2 + (3 2 )(r) = . (4.16) 

Proof: As in the proof of Proposition 14. 3\ we consider the Maxwell equation (|4.14|) . Since 
\4>'\ is bounded near the horizon according to condition (II), we conclude from (I) that the 
inhomogeneous term in (|4.14p must be integrable, 

/ P+ ^(« 2 + /? 2 ) < oo . (4.17) 
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Next we take the derivative of the function in (|4.16[) 
d 



dr 



(r-p)-2 (a 2 +P') 



S -(r- p)-^ (a 2 + (3 2 ) + (r - p)~i |-(a 2 + (i 2 ) , 



and substitute the bound (|4.10p . 

d 
dr 



((r -p)"f (a 2 + /3 2 )) 



(2j + 1 



i 

,2\ 2 



A~3 (r-p)~3 (a 2 + /? 2 ) . (4.18) 



Since s > 2, we have (r — p) 2 1 < (r — p) s , and thus (|4.17p implies that the first 
summand on the right side of (|4.18p is integrable. Using 



A -\ ( r _ p )-f ( a 2 + ^2) = A -l (a 2 + ^2) \ A \ ( r _ /f)) - 



0(A- 1 (« 2 + / 3 2 )) , 



we see that, according to ()4.17j) . the second summand in (|4.18p is also integrable. As a 
consequence, the function (r — p)~^ (a 2 + (3 2 ) has a limit at r = p. If this limit were 
non-zero, the integral (|4.17j) would diverge. We conclude that this limit must be zero. ■ 



Lemma 4.5 If s > 2, the function \<f>'\ has a finite, non-zero limit on the horizon; namely 

lim \(j>'\ = - lim T _1 > . (4.19) 

p<r^p p p<r^p 

Proof: To begin, we first show that 

lim (lu — e4>) T 2 (a 2 + f3 2 ) = . (4.20) 

p<r^p 

If the function \{u) — e<f))T\ is bounded, then (|4.20p is an immediate consequence of Lemma 
14.41 Thus we must only consider the case that — e<f>)T\ is unbounded near the horizon. 
The differential equation for AT 2 , (|4.12p . gives the estimate 



d Jm9 
r —AT 
dr 



> T 3 (a 2 + (3 2 ) ( 4(2j + 1) I (a; - e<t>)T\ - 2 (2j + 1)2 - 2(2j + 1) m 



According to assumption (I), the left side of this inequality is bounded near the horizon. 
Using that \(co — ecf))T\ becomes arbitrarily large near the horizon, we conclude that the 
function T 3 (a 2 + (3 2 ) \(u — e<p)T\ must be bounded. This implies (|4.2U|) . 

We now consider the 4-equation (|3.1ip . Since s > 2, the left side of (|3.1ip converges 
to zero for r \ p. Thus the right side of (|3.1ip must also go to zero in this limit, 

0= lim 1- A- 2(2 j + - ecf))T 2 (a 2 + I3 2 ) -r 2 AT 2 \<p'(r)\ 2 . 
p<r—+p 

Using (|4.2Up completes the proof. ■ 
We can now rule out the case s > 2; namely we have 
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Proposition 4.6 For s > 2, there are no solutions of the system /I3.y\) ~ €3.1ty) . 



Proof: According to Lemma 14.51 the function (u — e<p) has a Taylor expansion around 
the horizon with non-vanishing linear term, 

(ui — e(j))(r) = c + d (r — p) + o(r — p) with \d\ = — lim T _1 . 

p p<r->p 

Thus the coefficient {lo — e(j))T in the Dirac equation (|3.9[) . (|3.10p is monotone near the 
horizon and diverges. Using f|3.9|) and (|3.10p . this implies that the vector (a, f3) spins 
around the origin faster and faster as r approaches the horizon, which suggests that it 
cannot go to zero in this limit. In fact, [31 Lemma 5.1] yields that the spinors are bounded 
away from zero, 

liminf (a 2 + 3 2 ) > . 

p<r^p 

This contradicts Lemma 1-4.41 ■ 



Proof of Theorem \4-l\ According to Proposition 14.31 and Proposition 14.61 we must only 
consider the case s = 2; thus (14. 2\\ and (14.31) hold. Lemma 14.41 yields that the function 
(a 2 + f3 2 ) goes to zero on the horizon. Applying [3j Lemma 5.1], one sees that the function 
(uj — ecj))T cannot diverge monotonically near the horizon. On the other hand, Lemma [4.5 
shows that (a; — e<fi) has a Taylor expansion around the horizon with non-vanishing linear 
term. We conclude that (|4.4p holds, and that {oj — e<p)T has a finite limit on the horizon, 



lim (u — ed)T = A with |A| = — lim (i — p) A 

p<r^p y ^' p p<r^p y 



Exactly as in Section 5], one can rewrite the radial Dirac equations as ODEs in the 
variable 

u(r) = —r — p ln(r — p) 

and apply the stable manifold theorem [H Thm. 4.1] to conclude that a and (3 satisfy a 
power law near the horizon, (|4.5p and (|4.6p . Lemma 14.41 gives the bound k > ^. 

Now we substitute the expansions (|4.2p - (|4.6p into the Dirac equations (|3,9p . (|3.10p . 
This gives the conditions 

2j + l 

2p 

2j + 1 

VA)kA) = -(efoTo + m)a =F — O , 

2p 

which are equivalent to (|4.7p and (|4.8p . ■ 



Our main theorem gives strong restrictions for the behavior of black hole solutions 
near the event horizon. According to the condition k > g, the Dirac wave functions must 
decay so fast in the limit r \ p that they have no influence on the asymptotic form of 
both the metric and the electric field on the horizon. Namely, according to (|4.2p - (|4.4p . 
the metric and electric field must behave near the horizon like a vacuum solution, more 
precisely like the extreme Reissner-Nordstrom solution. The restriction to the extremal 
case means physically that the electric charge of the black hole must be so large that the 
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electric repulsion balances the gravitational attraction and prevents the Dirac particles 
from "falling into" the black hole. This is certainly not the physical situation which one 
can expect in the gravitational collapse of e.g. a star in the Universe. Nevertheless, extreme 
Reissner-Nordstrom black holes have zero temperature [6] and can thus be considered as 
the asymptotic states of black holes emitting Hawking radiation. For this reason, it is 
interesting to study if the asymptotic expansion of Theorem 14.11 leads to global solutions 
of the EDM equations. 

For an extreme Reissner-Nordstrom background field, it is proven in [3l Section 5] 
that the solutions of the Dirac equation satisfying (|4.5p . (|4.6p violate the normalization 
condition (|3.14p . Thus the question is if the influence of the spinors on the gravitational 
and electromagnetic fields can make the normalization integral finite. This is a hard 
analytic problem, because one must control the global behavior of the solution. However, 
we have done numerical investigations, taking the expansions in Theorem 14.11 as initial 
condition on the horizon and solving the equations for increasing r. It turns out that the 
solutions either develop singularities for finite r, or the spinors (a, (3) are not normalizable. 
Thus our numerics show that the expansions in Theorem l4.1l do not give normalizable black 
hole solutions of the EDM equations. 
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